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Abstract 



This article discusses the occurences of one-loop amplitudes within a next-to-next-to-leading 
order calculation. In an NNLO calculation the one-loop amplitude enters squared and one 
(^ ■ would therefore naively expect that the 0{z)- and o(£^) -terms of the one-loop amplitudes 

are required. I show that the calculation of these terms can be avoided if a method is known, 
which computes the O (£*') -terms of the finite remainder function of the two-loop amplitude. 
> ■ 

X 



1 Introduction 

Precision calculations in high-energy particle physics require the computation of higher order 
perturbative corrections. In the area of jet physics one needs in particular fully differential cal- 
culations which allow arbitrary infrared-safe experimental cuts. Several processes have been 
calculated in the past to next-to-next-to-leading order (NNLO) accuracy [1-29]. Most of these 
calculations have been done with an approach based on Feynman diagrams. 

At the level of next-to-leading order (NLO) calculations there has been a significant break- 
through for multi-particle amplitudes in recent years. This has been mainly due to unitarity or 
cut-based techniques [30-40]. In addition, purely numerical approaches based on subtractions 
for the computation of multi-parton one-loop amplitudes have been developed [41-44]. It is a 
common feature of these approaches that they no longer rely on Feynman diagrams but work 
on the level of amplitudes instead. It is a question of current interest if these methods can be 
extended to NNLO. This raises immediately a question related to the one-loop amplitudes. At 
NNLO the one-loop amplitude enters squared, and since the expansion in the dimensional regu- 
larisation parameter £ starts at order (—2) one would naively expect that up to order e^ the O (e)- 
and (£^)-terms of the one-loop amplitude are needed for an NNLO calculation. In fact, several 
groups (including the author) have calculated 0{e)- and o(£^) -terms of one-loop amplitudes in 
the past [45-53]. However, it is by no means obvious how the new approaches for one-loop am- 
plitudes based on unitarity or subtraction can be extended to include the higher-order terms in the 
£-expansion. For the unitarity approach the complication arises through the required generalisa- 
tion of the so-called rational terms beyond O (e^), while the numerical method with subtraction is 
based on evaluating finite integrals in four space-time dimensions, and thus is a priori insensitive 
to the additional (—2£) -dimensions. 

It is therefore appropriate to investigate first what is really needed for an NNLO calculation. 
This is the purpose of this article. In this paper I will trace every appearance of the one-loop 
amplitude inside an NNLO calculation. I will show that the O (£)- and (£^)-terms of the one- 
loop amplitude drop out from the final result provided a method is known, which computes the 
(£*^)-terms of the finite remainder function of the two-loop amplitude. Therefore what is actu- 
ally needed is just a method to compute the (£'^)-terms of the finite remainder function of the 
two-loop amplitude. In view of the new techniques for the computation of one-loop amplitudes, 
this is conceptually simpler than the computation of the two-loop amplitude up to order O (e^) 
and the computation of the one-loop amplitude up to order (£^) . 

This paper is organised as follows: In the next section the problem is exposed. Section 3 in- 
troduces the notation for this article and reviews the known results, which will be relevant in the 
following. The topics covered in this section are the ultraviolet renormalisation of loop ampli- 
tudes, the infrared pole structure of loop amplitudes, the behaviour of tree and loop amplitudes 
in singular regions of phase space as well as a short introduction into the subtraction method. 
Section 4 is the main section of this article. In this section I derive a theorem, which states that 
only the tree-level amplitudes, as well as the (£*^)-terms of the finite remainder functions jF ^^' 
and 5 ^^' are needed for an NNLO calculation. Finally, section 5 contains the conclusions and an 
outlook. In an appendix the changes for going from Catani's original definition of the insertion 
operators and remainder functions to a definition in a minimal scheme are discussed in detail. 



2 Exposition of the problem 

The essential point of the argument can already be explained for the case where all infrared 
singularities occur in the final state. Therefore I will discuss the case of electron-positron annihi- 
lation first and in detail. The case of hadronic initial states requires only minor modifications and 
is explained at the end of section 4. The master formula to calculate an observable at an collider 
with no initial-state hadrons (e.g. an electron-positron collider) is given by 
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^^) = ^^FTT 7T\ TT^L / ^4'nO„ (/?!,...,/?„, ^1,^2) l-^nP, (1) 

2i^(5)nspin(l)nspin(2) V-^ 

where q\ and ^2 are the momenta of the initial-state particles, 2K{s) = 2s is the flux factor 
and s = {q\ + qiY is the centre-of-mass energy squared. The factors l/nspin(l) and l/nspin(2) 
correspond to an averaging over the spins of the initial particles. d(^n is the invariant phase space 
measure for n final state particles and 0„ (;?i, ...,;?„, ^1,^2) is the observable, evaluated with a 
configuration depending on n final state partons and two initial state particles. The amplitudes 
^n are calculated perturbatively. The leading-order contribution to a n-jet observable is given by 

\^ |2 _ ^(0)*_^(0) .2) 

At NLO we have the virtual and the real contribution: 

I '^INLO — " " " " ' 

l^^,|2 _ (0) * (0) .3 

At NNLO we have the following contributions: 
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I -^"+2 1 NNLO ~ -^11+2 -^n+2- V^; 

Here J^n denotes an amplitude with n final-state partons and / loops. We see that in the NNLO 
contribution the one-loop amplitude with n final- state partons enters squared and that the one- 
loop amplitude with (n + l) final-state partons enters interfered with the corresponding tree am- 
plitude. Usually the loop amplitudes are calculated within dimensional regularisation. I denote 
by £ = (4 — D) /2 the regularisation parameter. The loop amplitudes have a Laurent series expan- 
sion in the parameter £. This series starts at order (—2/) for an Z-loop amplitude. In particular, 
for an one-loop amplitude one has 
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The square of the one-loop amplitude up to order (e*^) is then 
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One would therefore naively expect that the £-part Rn ' and the £ -part ^„ ' of the one-loop 
amplitude Rn are needed for an NNLO calculation. In this paper I will investigate if the calcu- 
lation of these quantities can be avoided. 

A related question concerns the all-order £-behaviour of the one-loop amplitude R^'^^^. 
NNLO -^„_j_i enters the contribution of the {n + 1) -particle final state 
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2i^(5)nspin(l)nspin(2) 

This expression contributes to an n-jet observable, whenever the (n + l) final state particles 
are classified by the experimental resolution criteria inside the observable O as an n-particle 
configuration. This happens when one of the particles becomes soft, or when a pair of particles 
becomes coUinear. The phase space integral is actually divergent, and the divergence results from 
the region where a massless particle is soft or where two massless particles are coUinear. This is 
the same situation as in an NLO calculation, where the real emission contribution is given as a 
phase space integral over (n + l) final state particles, weighted by the matrix element 
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In both case one usually applies the subtraction method [54-59] to handle this problem. Within 
the subtraction method one subtracts and adds back a suitable chosen approximation term. The 
approximation term must have in all singular limits the same behaviour in D dimensions as the 
approximated matrix element. In this way, the difference between the matrix element and the 
approximation term is integrable over the complete phase space. In addition the approxima- 
tion term should be simple enough, such that a phase space integration for one particle can be 
performed analytically. After this integration has been performed, the expression "lives" on an 
n-particle phase space and is combined with the other terms contributing to the integral over the 
n-particle phase space. 

In the beginning of the era of NNLO calculation it was feared, that one would need to know 
the all-order £-behaviour of ^„^i in order to find a subtraction term, which approximates this 
amplitude in D dimensions in the singular limits. Luckily, it was recognised soon that a one-loop 
amplitude factorises in the singular limits in a universal way of the following form [60-66] 
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The notation Sing o j?„ indicates that the factorisation occurs at the level of primitive amplitudes. 
Sing^*^) and Sing*^^-* are universal functions describing the singular limit (soft or coUinear). The 
universality of the singular functions Sing'*^-* and Sing'^-* implies that they do not depend on the 
complete set of momenta {pi, ...,/?„}, but only on two or three external momenta. In particular 
they do not depend on the number n. In order to find the approximation terms, only the D- 
dimensional behaviour of the singular functions Sing*^*^) and Sing*^^) need to be known, not the 

D-dimensional behaviour of the one-loop amplitude ^„_^i ■ Therefore, for an NNLO calculation 
one needs to know only the first three terms of the £-expansion of the one-loop amplitude ^„^i : 

4':. = ^^!^" + l^!^"+^!^ + o(e). (10) 

This fact is well established [60-66]. 

One observes that in the factorisation formula eq. (9) the one-loop amplitude Rn with n 
particles in the final state appears. In addition, there is a third place where the one-loop amplitude 
Rn appears: The two-loop amplitude Jl„ can be written as a part containing all the poles and 
a finite remainder. There is a universal formula due to Catani [67], which describes the pole part. 
In this formula the one-loop amplitude J^„ enters. 

In the following I will show that in the combination of the three occurrences of the one-loop 
amplitude Jin the 0(e)- and o(£^) -terms drop out. 

3 Notation and review of known results 

3.1 Ultraviolet renormalisation of loop amplitudes 

The loop amplitudes have explicit divergences. It is common practise to regulate these di- 
vergences by dimensional regularisation. Within dimensional regularisation one continues the 
space-time dimension from four to D = 4 — 2£. The origin of these divergences are either related 
to ultraviolet or to infrared singularities. The ultraviolet divergences are removed by redefining 
the parameters of the theory. In massless QCD it is sufficient to renormalise the strong coupling. 
In the MS scheme the relation between the bare coupling ao and the renormalised coupling 
(Xs{iJ^) evaluated at the renormalisation scale /j^ reads: 



ao = a,SgV^'' 
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where 
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is the typical phase-space volume factor in D = 4 — 2£ dimensions, Je is Euler's constant, and (3o 
and Pi are the first two coefficients of the QCD (3-function: 

Po = Y^A - ^TRNf, pi = y Ci - ^CAlRNf - ACpTRNf, (13) 



with the colour factors 
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Let the expansion in the strong coupling of the unrenormalised amplitude be 
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Then, the renormalised amplitude can be expressed as 
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The relations between the renormalised and the bare amplitudes are given by 
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In this paper I work with renormalised amplitudes and I will drop the subscript "ren" in the 
following. 



3.2 Infrared structure of loop amplitudes 

The infrared pole structure of loop amplitudes in the dimensional regularisation parameter £ is 
well understood [67-78] and can be stated explicitly for one- and two-loop amplitudes. If we 
regard a loop amplitude as a vector in colour space, then the infrared poles of the loop amplitude 
can be expressed through an operator acting on this vector in colour space. It is therefore conve- 
nient to introduce the colour charge operators T,. The action of a colour charge operator T; for a 
quark, gluon and antiquark in the final state is given by 



quark: ^* (...^,...) (t;.^) ^ (...^,...) , 
gluon: ^* (.../...)(//-*) ^ (.../. 
antiquark : SI* {...qi...) {-Tf^ s^ (•••<?/•••) • 



(18) 



The corresponding formulae for colour charge operators for a quark, gluon or antiquark in the 
initial state are 



quark: J?* (...^,...) (-T,^) ^ (...<?,-...) 
gluon: ^* (.../...)(//-*) ^ (.../... 
antiquark: Si* {...qi...){T['j) SA [...qj...) . 



(19) 



In the amplitude an incoming quark is denoted as an outgoing antiquark and vice versa. We start 
with the one-loop amplitude, which can be written as 

^^ = I«^r + iF„«. (20) 

Here l(^) contains all infrared double and single poles in l/£ and ir„ ' is a finite remainder. The 
e-expansion of ir„ starts at order t^: 

ir« = iFi^'°)+E5i^'^) + eVi^'^) + o(e3). (21) 

At two-loops, the corresponding formula reads: 

^P = i(2)^r + i(i)^ii) + :ri^). (22) 

Again, the £-expansion of the remainder function starts at order £^: 

jP = iFi^'°)+85i^'^)+eVi^'^) + o(e3). (23) 

One observes that in eq. (22) the one-loop amplitude Rn occurs in combination with the one- 
loop insertion operator I^^). The one-loop insertion operator I^^^ is given in the massless case 
by 



where 



I'iie) = T^l+yi, (25) 



and the coefficients T? and y, are 



t2 = T? = Cf, t2=Ca, y, = y^ = ^C^, Yg = ^. (26) 

In general, the colour operators T/Tj give rise to colour correlations. The corresponding for- 
mula for the two-loop insertion operator I^^^ is known, but not needed in the present article. The 
generalisation to the massive case is also known [67-78]. It should be noted that eq. (24) corre- 
sponds to the original definition of the insertion operators I^^^ and I^^^ due to Catani. Within this 
definition, the insertion operators contain apart from the pole terms also terms of order £* with 
k>0. Alternatively, it is possible to define the insertion operators in such a way that they contain 
only the pole terms. This is discussed in detail in appendix A. It should be kept in mind that 
a different definition for the insertion operators will also change the finite remainder functions 

:r(i)andir(2)_ 



3.3 Singular behaviour in phase space 

QCD amplitudes become singular in phase space, when the momenta of one or more external 
particles become degenerate. In perturbative calculations this phenomen occurs first in next- 
to-leading order calculations. Singularities occur if either two particles become coUinear or if 
one particle becomes soft. The factorisation properties are most conveniently discussed through 
decomposing QCD amplitudes into primitive amplitudes. We may write 



^i" = Ici>«. (27) 



where the coefficients C^ ■ carry all the colour information. The quantities A^ are called primi- 
tive amplitudes. Primitive amplitudes have a fixed cyclic ordering of the QCD partons, a definite 
routing of the external fermion lines through the diagram and a definite particle content circulat- 
ing in the loop. In the following I will discuss the factorisation properties of a single primitive 
amplitude and I will therefore drop the subscript j. For tree-level amplitudes one has in the soft 
limit the factorisation 

^^^ ^n+li-^Pi^Pj^Pk^-) = ^^^ iPi^Pj^Pk)An i-^PuPk^-): (28) 

Pj soft 

where the eikonal factor is given by 

Eikf(p„p„n) = ^J1<£A^}£^<£A. (29) 

The square of the eikonal factor is given by 

Y^Eikf\pi,pj,pkfEikf\pi,pj,pk) = 4^^. (30) 

Soft singularities lead to colour correlations among the primitive amplitudes A„ '. 

In the coUinear limit one has the factorisation 

limAi'Ji(...,p,p;,...) = Y^Sp\itf\p,,pj)A^^\...,p,...). (31) 

PiWPj A, 

where pi and pj are the momenta of two adjacent legs and the sum is over all polarisations. 
Squaring the splitting amplitudes one obtains 

<uark(^. ^0 = I "(p. ^) Splitf (P„ Pj) * Splitf {pu Pj) U{p, I') , 

^Suon(^A') = l^e'ip^r Spmf\pi,p/Sp\itf\pi,pj)e'{p^'). (32) 



P^^' {X, X') is a tensor in spin space. CoUinear singularities lead therefore to correlations in spin 
space. The spin- averaged splitting functions are obtained by 

,(0) \ _ 1vp(o) 



V3,quark) ~ 9 1^^3,quark(^'^) 



The corresponding formulae for the factorisation of one-loop amplitudes in soft and coUinear 
limits are [60-66] 

J soft ^ 

Eikf\pi,pj,pk)Ai^\...,Pi,Pk,...)+'Eik^\pi,Pj,Pk)A''n\...,Pi,Pk,...), 

lirnAl^l^{...,Pi,Pj,...) = 

PiWpj 

£ (Sp\itf\pi,pj)A^^\...,p,...) + Sp\it^i\p,,Pj)AT{...,p,...)). (34) 

The one-loop singular functions Eikg and Splitg are known [60-66]. In this article we are 
concerned with the part proportional to Aa . Here the tree-level singular functions Eikg and 
Splitj occur. Therefore this part is very similar to the NLO case. 

3.4 The subtraction method 

The individual contributions to {o)nlo and {o)nnlo are in general infrared divergent, only the 
sum is finite. However, these contributions live on different phase spaces, which prevents a naive 
Monte Carlo integration approach. To render the individual contributions finite, one adds and 
subtracts suitable chosen terms. At NLO we have [54-59] 

{0)nLO = j (^0n+ldO^^l^-0nOda^n'^^^+ J {OndO^n^ + OnOda^n'^^y (35) 

The notation 0„ o ddn is a reminder, that in general the approximation is a sum of terms 

OnOdOh = ^Ondan (36) 

and the mapping used to relate the n + \ parton configuration to a n parton configuration differs 
in general for each summand. 

In a similar way, the NNLO contribution is written as [79-85] 

( o)nnlo = / yOn+2 dofl2 - 0„+i o da^^^} - 0„ o Jai°'^^ 
+ / {on+\ do[^l^ + On+1 o dal^_l^^ - o„ o dai^ 
+ I i^On d(5^n^ + Onodaf'^^ + o„o Jai^'^^j . (37) 
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(iajj^j is the NLO subtraction term for {n+ l)-parton configurations, dai ' and dali ' are 
generic NNLO subtraction terms. These terms can be decomposed further into the following 
form 

dai!'^^ = Ja|,°°P + Jar'^"''-rfa^i°^°''-J<°" + Jar^"='^. (38) 

The terms Ja™^ ^ and d(Xn°^ are of relevance here, da^™^ ^ is the NLO subtraction term. The 
term <ia^°°P approximates the singular behaviour of one-loop amplitudes in the soft and coUinear 
limits. This term can be written as the sum of two contributions 

Ja'°°P = Ja'°°P'" + Ja'°°P'^ (39) 

such that (ia'°°P'" contains the tree-level approximation functions and one-loop amplitudes, while 
(ia'°°P'^ contains the one-loop approximation functions and tree-level amplitudes. In formulae 
we have 

lim Ja'°°P''^ = Sing(°)oJai^\ 

n+l— >n 

lim Ja'°°P'^ = Sing^^^odol^l (40) 

n+l— >-w 

The one-loop amplitudes A„ appear in (ia^°°P''^ and we will focus on this term in the sequel. 

We need the integral over the unresolved phase space of Ja™^ ^ and 6?a„°°P'". For the NLO 
subtraction term we have 



/"""'<'" = ik) n^Mn,U2) I '''■"'■ ^''''i''^'''-'^^'^'' <^" 



where the operator l^J^^ contains all the poles in the dimensional regularisation parameter £. 
^reai ^^ independent of the subtraction scheme. F^^^j is a finite remainder and depends on the 
subtraction scheme. I^^^j is given in massless QCD by 

'-1 - 2;rr(l-8)-V,^ ' ^U^^T^eJI^ ^^ ) ' ^^'^ 

Note that 1^^!^^ is very similar to I*^^'. The insertion operator 1^^!^^ is obtained from I^^^ by multi- 
plying with a factor (—2) and by the substitution 

-IpiPj -^ \^PiPj\- (43) 

For the integral over the unresolved phase space of dCLn°^'"' one has 
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4 Cancellation of the higher-order terms 

We are now in a position to put all pieces together. It is instructive to study the NLO case first. 
For the virtual part we have 

2Kc^^^\l^^ = 2Re^r*l(i)^r + 2Re^rVi^) (45) 

The integration of the subtraction terms for the real part adds a term 

In the sum we have 

Re 4"' ' (2I<" +lS) ^i°' + 2 Re 4" V„"' + ^f 'pW ^f . (47) 

The higher-order terms of the one-loop amplitude are contained in the remainder function ir„ . 
This function gets interfered with tree amplitude Jl„ , therefore only the first term J„ ' of the 
e-expansion is needed. It is worth noting that the combination 21^^) +Ireai ^^ ^^ general not free 
of poles: A single pole remains, whose coefficient is purely imaginary: 

^'^+^ = f|llT,T^-0(-2w) + o(e°)- (48) 

This pole drops out by taking the real part, therefore eq. (47) is free of poles. 

Let us now move to NNLO. We consider all terms, where the one-loop amplitude Rn occurs. 
First of all, we have the square of the one-loop amplitude J^n . We can write this term as 

^'^ = ^r*I«*lO)^r+2Re^r*lO)V« + 5i^)Vi^). (49) 

Secondly, from the integration of <ia„°°^'" we get 

2 Re :, <" • (l« + fS) ^<- = 2 Re ^i»' 'iW I<>'^i°' + 2 Re ^i«' * I« :r,."' 

+2 Re ;,<» > •Fi:],I<');,i°' + 2 Re ^^f ' F^, /„<". (50) 
Thirdly, the contribution from the two-loop amplitude can be written as 

2Re^r*^P = 2Re^r*l(i)l(i)^r+2Re;ir*l(%i^) 

+2Re^r*l(2)^r+2Re^rVi^). (51) 

Adding these three pieces together we find 

2 Re ^r ' ^ + ^ * ^ + 2 Re ^ * (ij^, + F« ) ^^^ = 

Re ^f * (21(2) +21(1)1(1) ^ Id) * Id) +21(^,1(1) +2f(2,i(i)) ^f 

+2Re^r*(l(i) + l(i)*+i;2,)5i^) 

+2 Re ^r * ^i^) + ^i^) * ^P + 2 Re ^f * F;2,:r„«. (52) 

11 



The terms in the first line of the r.h.s of eq. (52) involve only the Bom amplitude J?„ and the 

insertion operators, but do not involve the one-loop finite remainder function ir„ . Let us then 

(1) * (1) 
consider the terms in the third line of the r.h.s of eq. (52). It is clear that for the terms ir„ J„ 

and 2 Re J^n ^y^li^n only the £*^-terms of J„ are needed, since J„ is multiplied by terms 
which do not have any poles in £. It remains to examine the terms in the second line of eq. (52). 
These are the most critical ones. They are given by 

2Re^f*(l(i) + lWVlS):rii^ (53) 

Although the individual insertion operators start at o(£ ^), it turns out that in the combination 
of the insertion operators all poles cancel 

l(i) + l(i)*+l£ii = 0(e°), (54) 

and it follows, that again only the £^-term iT,) ' of the one-loop remainder function is needed. 
We illustrate the cancellation of the poles for the massless case. Here we have the explicit 
expression 



i(i) + i(i)Vi!! 



real 



a, e^Y£ n ji 1 



27rr(l-£)^t^. '\e 



T2£ 



Re / ^P'Pj\ ' (\^PiPj\ 



1-^ ) I A*^ 



(55) 



The expression in the square bracket is of order £?, therefore the complete expression starts at 
order £^. 

For completeness let us also discuss which terms are needed for the finite remainder func- 
tions y^n and iF„_^{ . The two-loop remainder function jfn occurs only in the interference term 
with the Bom amplitude Jl„ , therefore it follows that only the £*^-term ^„ ' ' of the two-loop 
remainder function is needed. The one-loop remainder function J^^^ with in-\-\) partons oc- 
curs only in the interference term with the Bom amplitude ^„_|_\ and in combination with the 
subtraction term Ja'°°P. This contribution is evaluated in four dimensions, and therefore only 

the £*^-term 5„_|_\ of the one-loop remainder function is needed. 
We may summarise the situation as follows: 

Theorem: For an NNLO calculation it is sufficient to know the tree-level amplitudes ^„ , -^„_|_\ 

and ^„_/2' ^^ 0(£^)-terms of the one-loop finite remainder function jF„ and O^I^Ji, as well as 

the 0(£'^)-terms of the two-loop finite remainder function iF„ . The 0(£)- and o(£^) -terms of 

the one-loop finite remainder function Jn drop out from the final result and are therefore not 
needed. 

This is the main result of this paper. A few remarks are in order: 
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1. Although not explicitly shown in this paper, the statement above can be sharpened: Only 
the four-dimensional results (or equivalently only the o(£'^) -terms) of the tree-level amplitudes 

Rn , ^n+i ^^'^ ^n+2 ^^^ required. This follows from the Kinoshita-Lee-Nauenberg theorem, 
which ensures that all poles cancel in the end. Therefore only finite expressions multiply the 
tree-level amplitudes and it is sufficient to know them in four dimensions. 



2. It should be pointed out that the statement above does not say, that the one- and two-loop 
amplitudes up to 0{£P) are sufficient. What needs to be known are the o(£'')-terms of the finite 



(2) 



remainder function 5„ , and not the O (e) -terms of the two-loop amplitude J?„ 



(2) 



3. The statement applies also to processes with partons in the initial state. This can be seen as 
follows: For processes with initial-state partons, the integration over the unresolved phase-space 
corresponding to eq. (41) and eq. (44) will give additional poles proportional to the Altarelli- 
Parisi splitting functions. These poles are canceled by a coUinear counterterm arising from the 
factorisation into a hard scattering process and parton distribution functions. The coUinear coun- 
terterm is given by 



do' 



1 1 

= ^ / dxi / dx2 

flJlQ 

1 1 



£ f d%0nff,{xi)ff^{x2) |^„P. (56) 



2K{s) nspin(l)«spin(2)ncolour(l)«colour(2) ^ 

The first sum is over all partons species /i and /2 inside the initial-state hadrons. //i(^i) and 
//2 (^2) are the parton distribution functions. The quantity 5 is the partonic centre-of-mass energy. 
The averaging is now also over the colour degrees of freedom of the initial state partons, given 
by ncoiour(l) and «coiour(2). Expanded to next-to-leading order, the coUinear counterterm yields 
a contribution proportional to 



Us 

2k^ 



5e / dzi / dz2\ 5(1 -Z2) 



U4 



+5(1 -zi; 



U4 



A* 



P^'\Z2)+K^'>{Z2 



A* 



(1) 



P^'Hzi)+K^'Hzi 



,(0) 



Jl„ ' {pi,...,Pn,ziqi,Z2q2, 



(57) 



P^^' are NLO Altrarelli-Parisi splitting functions. The functions K^^ are finite and define the 
factorisation scheme. In the MS-scheme one has ^*^^^ = 0. At NNLO the part involving the 
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one-loop amplitude J^n is very similar: 

1 1 



■|5e/jzi/jZ2 5(1-Z2) -^(^) P^'Hzi)+K^'Hzi: 



+5(l-zi) 



-i{§yp''Hz2)+K^'Hz2) 



z is£ yifi -H-Yi , 



(58) 



where it should be understood that the amplitudes are again evaluated with the rescaled momenta 
z\q\ and ziqi for the initial-state partons. The combination with the additional poles from the 
integration over the unresolved phase-space corresponding to eq. (41) and eq. (41) yields in the 
NLO case the finite contribution 



i 1 

y^zi I Jz2 {5(1 -Z2) [P(') (zi) +k(i) (zi)] +5(1 -zi) [p(i) (z2) +k(i) (z2)] } 



-^n 



(0) 



.(59) 



The colour-charge operators P^^-' and K(^) are finite. Their precise definition can be found in 
ref. [56]. In the same way one obtains at NNLO for the part involving the one- loop amplitude 



^n 



(1). 



1 1 

dz\ \ dzi 





(60) 



{5(1 -Z2) [P('^ (zi) +k(i) (zi)] +5(1 -zi) [p(i) (Z2) +k(i) (z2)] } 2Re Jli°^*^i'^. 

Since P*^^) and K(^) are finite, the 0{€)- and o(£^)-terms of the one-loop amplitude are not rel- 
evant for this contribution. This completes the proof of the theorem for the case of initial-state 
partons. 

At the end I would like to discuss a possible application: Turning eq. (20) and eq. (22) around, 
we have 



^(1) _ a(i)_i(i)^(o) 



7n 



(2) 



a(2)_i(l).(l)_i(2) (0) 



(61) 



The remainder functions Jn and J„ are finite quantities. The (e^) -terms can be calculated 
numerically in four dimensions as follows: One first introduces an integral representation for the 
insertion operators l(^) and I^^^ in D dimensions. The right-hand side can then be viewed as a 
single integral over 



(27rp 
14 



(62) 



for fn and as a single integral over 

(2) 

for ir„ . In the next step one adds an approximation term, which approximates the integrand 
locally in all singular limits and which integrates to zero. In the final step the limit D — t- 4 can 
be taken and the integrals can be evaluated in four dimensions. In particular the integrals can 
be evaluated numerically by Monte Carlo methods. At NLO this approach has been studied 
in [41-44]. The present article shows that this method can be generalised to NNLO, bypassing 
the need to calculate the (e)- and O (e^) -terms of the one-loop amplitudes. 

5 Conclusions 

In this article I traced the appearance of one-loop amplitudes within an NNLO calculation. I 
investigated whether the calculation of the 0{e)- and C'(£^)-terms of the one-loop amplitudes 
can be avoided. I showed that the calculation of these terms can be avoided, if a method is 
known which allows the computation of the O (£*^) -terms of the finite remainder functions for the 
two-loop amplitude and the one-loop amplitude. This opens the door for a numerical calculation 
at NNLO with the help of the subtraction method. The result of this paper will also be useful for 
NNLO calculations within the unitarity or cut-based method. 
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A Scheme dependence of the insertion operators 

In the main text of this article I have used the original definition of the one-loop and two-loop 
insertion operators 1*^^^ and I^^^ due to Catani. With this definition the insertion operators do 
contain apart from the pole terms also terms of order £^ with A: > 0. It should be mentioned that 
these insertion operators can also be defined in a minimal scheme, such that in this scheme they 
contain exactly the pole terms and nothing else [70,71]. I will denote the insertion operators in 
the minimal scheme by Z^ '> and Z^ K The relation between the two schemes is given by 



Z(i) = r[iW 
Z(2) = 7?(l(2)+l(i)7?(l(i))), (64) 

where R denotes the projection onto the pole part: 

R f I 'jA = t 'J^'- (65) 
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Changing the insertion operators to the minimal scheme will also change the finite remainder 
functions. One has 



.(1) _ zW^^'^)+crW 

Jin —Li -^n ^-r„m\ 



n, minimal' 



,m ^ ^z<«-z(»z<>));,f +z<');,<'' + :r,S.„„., (66) 



The result of this paper is not affected by this change of schemes. Therefore the knowledge of 

the (£°) -terms of the 
an NNLO calculation. 



the (£^) -terms of the finite remainder functions ^„ ^^inimai ^^'^ ^n minimal ^^ equally sufficient for 
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